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Abstract. By considering the least-squares fits of the 
double-mode Cepheid light curves discussed in Paper I 
wc defined their properties by their Fourier parameters 
and generalized phase differences G, j . When plotting the 
latter quantities as a function of the respective order, the 
second order terms are confined in the region just below 
3/27r; the third order terms have ir/2 < Gi.j < tt, the 
fourth order ones cluster around 2tt (or 0), the fifth order 
ones seem to have ir < Gj j < 3/2tt. The mean G^j values 
are also regularly spaced. The progression of the G,-j val- 
ues as a function of the period was investigated and the 
signature of a possible resonance near 6.0 d was found. 
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1. Introduction 

Pardo & Poretti (1996, hereinafter Paper I) made a fre- 
quency analysis of the available photometry of galactic 
double-mode Cepheids (DMCs) and obtained a very re- 
liable set of Fourier parameters for each star. However, a 
fully comprehensive and synthetic observational descrip- 
tion of the DMC light curves is still lacking. In this work, 
directly originating from the previous one, we try to deter- 
mine the common characteristics of the DMC light curves, 
searching for the boundary values of the Fourier parame- 
ters and regularities in their progression. 



2. The generalized phase differences 

In Paper I we fitted the V magnitudes by means of the 
formula 



V(t) = V + }_^A Z cos[2ir f z (t - T ) 



(1) 



where f z is the generic frequency, which can be an in- 
dependent frequency (/iand f2), a harmonic or a cross 
coupling term. Our previous analysis demonstrated that 
each component in the DMC light curves can be defined as 
a combination of two basic frequencies /iand fi\ by defin- 
ing z = we have f z — fij=ifi+jf2- Some examples: 
for (i,j)=2,0 we have the harmonic 2/i; for (i,j)=l,l the 
fi + pterin; for (i,j)=-l,l the f 2 - Jitcrm; for (i,j)=3,2 
the 3/i-/2and so on. 

In order to define the properties of the Fourier param- 
eters of the DMC light curves it is very useful to recall to 
mind the generalized phase differences introduced by An- 
tonello (1994b A&A 291, 820), here noted as G id . They 
are a linear combination of the phases of each term fo j 
and of the phases <&i and $2 of the independent frequen- 
cies /iand /2(^i and ^2 in Antonello's notation). Their 
expression is given by 



G 



^,3 



J h3 



(2) 



The numerical application to the TU Cas fit provides 
some examples (the integer k values have to be selected 
so that G id G [0,27r]): 



Gi, 



h,i - $1 - $2 



2fc7T 



4.13 - 4.31- 2.19 + 2tt = 3.91 



G_l I = 0_i 1 + $1 - $2 + 2fc7T = 



= 1.92 + 4.31 - 2.19 = 4.04 
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G 3 , 2 = 03,2 - 3$i - 2$ 2 + 2&7T = 



2 
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Table 1. Generalized phase differences for all the galactic DMCs; their measure units are in rad 10 -2 . The amplitude ratios 
R21 for the /iand ^frequencies are also reported. Except for the unique 10/20 pulsator CO Aur, the stars are listed in order 
of increasing period 



TU Cas U TrA VX Pup AS Cas AP Vel BK Cen UZ Cen 



R-2l(/l) 


346±4 


323±4 


184±6 


R-2l(/ 2 ) 


123±9 


99±20 


118±7 




415±3 


415±3 


414±8 


/1 + /2 


391±4 


404±3 


414±5 


/2-/1 


404±6 


374±3 


450±9 


2/2 


433±10 


439±10 


447±13 


3/1 


207±6 


208±6 


231±33 


2/1 +h 


189±6 


246±5 


206±13 


A +2/2 


178±11 


181±9 


228±17 


2/2-/1 


165±21 




238±25 


4/i 


604±13 


588±13 




3/i +A 


608±9 


623±8 


668±34 


2/1+2/2 


586±13 







3/1-/2 

3/1+2/2 364+21 

4/1+/2 386+15 407+15 



305+10 


283+4 


263+8 


326+11 


139+15 


117+15 


102+19 


122+37 


415+5 


417+4 


423+5 


420+6 


412+4 


396+5 


427+7 


416+11 


409+7 


441+6 


445+11 


393+20 


434+11 


445+13 


488+23 


449+39 


219+15 


231+9 


208+11 


216+10 


210+8 


222+9 


214+9 


226+14 


243+12 


202+17 


230+23 




219+18 














642+26 


629+19 


588+27 


621+24 


634+21 



687+52 



Y Car AX Vel GZ Car BQ Ser EW Set V367 Set CO Aur 



R2i(A) 


298+12 


104+19 


156+14 


R2i(A) 


103+26 


77+7 


46+23 


2/1 


418+6 


417+9 


433+10 


A + A 


427+7 


436+5 


453+10 


A -A 


458+7 


466+9 


428+14 


2/2 


433+13 


456+10 


470+49 


3/i 


226+13 






2/1 +A 


224+17 


283+18 




A +2/2 


227+34 






2/2-/1 






165+25 


4/i 








3A +A 


638+30 







175+6 


164+12 


153+6 


179+6 


45+9 


24+9 


120+9 




421+5 


440+7 


450+7 


409+8 


450+4 


469+8 


402+9 


452+24 


456+6 


469+10 


499+11 


473+28 


523+21 


300+34 


362+12 




230+19 


255+53 


218+36 


183+29 


234+13 


234+33 







= 2.11 - 3 • 4.31 - 2 • 2.19 + 6tt = 3.64 

Table 1 lists the Gij values for each star's Aj's, as 
calculated on the basis of the phase values listed in Tab. 
3, 4 and 5 of Paper I. For each value the maximum error is 
also reported, as obtained from the formal errors reported 
in Paper I. We remind the reader that the real significance 
of a component was established star by star and hence we 
are dealing here with very reliable values. Table 1 supplies 
the first synthetic description of the DMC light curves, as 
obtained from the available photometry. 

2.1. Separations between different order terms 

It is quite interesting to plot the Gij values against their 
fit order. Light curves of DMCs are often quoted as an ex- 



ample of erratic behaviour and cycle-to-cycle variations, 
both in amplitude and in phase. In Paper I we already 
proved that these light curve seem to be much more sta- 
ble than reported and that a frequency locked fit yields 
a satisfactory representation. Only in the cases of U TrA 
and EW Set we found some slight evidence of frequency 
or amplitude variations. 

The suspicion that the DMC light curves have a pre- 
dictable behaviour is confirmed by the natural upper and 
lower limits that can be easily observed in Fig. 1. The 
second order terms are confined in the region just below 
3/27r; the third order terms have ir/2 < Gij < n, the 
fourth order ones cluster around 27r (or 0), the fifth order 
ones seem to have ir < Gij < 3/2n. 
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The mean Gij values are 4.30±0.34 rad for the second 
order (i.e. \ i | + | j =2), 2.20±0.23 rad for the third one, 
6.24±0.31 for the fourth one, 3.85±0.21 for the fifth one. 
These mean values are roughly equispaced, with a slight 
tendency to increase: indeed, the differences between the 
mean Gij of adjacent orders are 2.10, 2.24, 2.39 rad, re- 
spectively. The latter result and the boundary values es- 
tablished above yield an experimental confirmation of the 
conjectures first expressed by Antonello (1994b) about the 
extension to DMCs of the rule of uniformity of phase dif- 
ferences in monoperiodic Cepheids. However, the observed 
separation (~2.2 rad) is a bit larger than expected (n/2) 
in the case of adiabatic pulsations in a one-zone model. 
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Fig. 1. A well defined separation is obtained when plotting 
the generalized phase values Gij in function of the respective 
order 



2.2. The Gij progressions 

The second order G itj values (2/i, /i + / 2 , fi - /i, 2/ 2 ) 
range from 3.00 to 5.23 rad; it was expected to see a lit- 
tle spread of the Go j2 values owing to the resonance at 
3.0 d. Indeed the two extrema are just related to the 
2/ 2 components of the BQ Ser (the DMC approaching 
resonance from the shorter periods) and EW Set (the 
DMC approaching resonance from the longer periods) 
light curves. Antonello (1994a) reported another possible 
resonance between the third overtone and the fi + / 2 term 
near 6.5 d. Later, Antonello (1994b) discussed a prelimi- 
nary progression of the Gi i values and stressed the im- 
portance of verifying the position of the points related 
to V367 Set and EW Set light curves. This can now be 
carefully done as Fig. 2. In the lower panel the last point 



Fig. 2. The cross-coupling terms /i + /2and /2 — /iare present 
in all the DMC light curves. While the G_i,i progression does 
not show any particularity, the Gi,i progression seems to be 
interrupted at P ~6.0 d 



(4.02±0.09 rad, V367 Set, P=6.293 d) is clearly out of 
the progression followed by the other points (the last is at 
4.69±0.08 rad, EW Set, P= 5.823 d). Moreover, the pro- 
gressive weakening of the amplitude of the /i + / 2 term is 
clearly visible in Fig. 3 (lower panel); the same trend can 
be evidenced by considering different types of normalized 
amplitudes and mode energies and it must be considered 
as a well established fact. On the other hand, if we look at 
the / 2 — /iterm we observed a smoother behaviour consid- 
ering both the G_i4 progression (Fig. 2, upper panel) and 
its amplitude (Fig. 3, upper panel). All these facts strongly 
supporting the action of a resonance effect involving the 
fi + / 2 tcrm. 

The third order progressions (G 2j i, Gi i2 , G_i, 2 ) do not 
show any particular feature except a slight tendency of 
increasing values for longer periods. The fourth order term 
G31 progression mimics, at a different level, the G3,o one. 



4 
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Fig. 4. The Fourier parameters fai and R21 provide a powerful discrimination between F and IO pulsators. Dots: single-mode 
Classical Cepheids. Triangles: s-Cepheids. Filled dots: Fundamental radial mode of DMCs. Filled triangles: IO radial mode of 
DMCs 



3. Double mode, classical and s Cepheids 

Figure 3 represents the actual scenario of single- and 
double-mode Cepheids with P <8 d in our galaxy. The 
4>2\-P plane provides a quantitative representation of 
the Hertzsprung progression, formed by the Classical 
Cepheids (open circles); the 02, 1 points obtained from the 
F mode of DMCs (filled circles) are perfectly superim- 



posed. The "Z" progression is defined by the single-mode 
Cepheids which do not follow the Hertzsprung progres- 
sion; in our opinion, these stars are to be considered the 
true s-Cepheids and their IO nature is now definitely es- 
tablished by the perfect supcrimposition of the fei values 
obtained from the 10 mode of DMCs. 
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ating only when considering the period, but they arc in 
the range of those of F/IO pulsators (see Fig. 2 for the 
fi + /2and fi — /i cases). In Paper I we already discussed 
its perfectly sine-shaped 20 light curve. 



4.0 

Period 



Fig. 3. The amplitude of the cross-coupling terms 
fi + /2 decreases towards longer periods and it reaches its mini- 
mum at P~6.0 d (lower panel); the amplitude of the fi — /iterm 
is small over the whole range of period values 



There is again a some underlying confusion about the 
nomenclature of these variables. The general assumption 
is that "s-Cepheids have a quasi-sinusoidal light curve", 
but to prove it it is necessary to perform a quantita- 
tive analysis: in such a case it is easy to establish the 
quasi-sinusoidal shape. Other criteria (first look evalua- 
tion, full amplitude values, ....) are quite arbitrary and 
not as useful as the amplitude ratios i? 2 i for /i(i.e. the 
ratio between the amplitude of /iand its harmonic 2/i) 
and /2 (between the amplitude of /2and 2/ 2 ). As the lower 
panel of Fig. 3 shows, the stars forming the "Z" sequence 
also show a small i?2i value and hence the light curve 
deviates very slightly from a sinewave shape. However, it 
should be noted that the R21 values for the F-mode of 
AX Vel (particularly), GZ Car, BQ Ser are smaller than 
the expected ones; also the R21 value of the .F-mode of 
VX Pup is small and it is located in the gap between the 
two sequences, very close to that of the single-mode AV 
Cir. Hence, the Fourier parameters have to be considered 
globally to perform a reliable mode identification. 

The Gij values of the fit of the CO Aur data were not 
used in the previous discussion because they strongly de- 
viate from the progressions described by the others; this 
is undoubtedly due to the 10/20 pulsating nature of this 
star. It should be noted that the dj values are devi- 



4. Conclusions 

In the previous sections we supplied a quantitative de- 
scription of the DMC light curves by means of their 
Fourier parameters and their progressions as a function 
of the period. We detected and determined very precise 
ranges of generalized phase differences in function of the 
order fit, also establishing the uniformity of the Fourier 
parameter distribution. In such a context, these unifor- 
mities make the DMC light curves very similar to the 
single-mode Cepheid ones. Hence, the physical properties 
of the pulsation should be the same and, at least to a first 
approximation, the adiabatic approximation provides a 
satisfactory understanding of the phenomenon (Antoncllo 
1994b) 

It will be very interesting now to verify the properties 
of the generalized phase differences in the more numerous 
sample constituted by the DMCs of the Large Magellanic 
Clouds (Alcock et al. 1996, Welch ct al. 1996); to do this, 
we consider it an essential step to apply the frequency 
analysis described in Paper I in order to avoid any spurious 
or unreliable result. 

Acknowledgements. The authors wish to thank E. Antonello 
for useful discussions and J. Vialle for the improvement of the 
English form of the manuscript. 

References 

Alcock C, et al. (MACHO project Team), 1996, AJ, 109, 1653 

Antoncllo E., 1994a, A&A 282, 835 

Antoncllo E., 1994b, A&A 291, 820 

Pardo I., Poretti E., 1996, A&A submitted (Paper I) 

Welch et al., 12 th IAP Colloquium, Paris, 1996 July, in press 



This article was processed by the author using Springer- Verlag 
UT E X A&A style file L-AA version 3. 



